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Abstract 

Recently, a new mechanism to generate a naturally small electroweak scale has 
been proposed. It exploits the coupling of the Higgs to an axion-like field and a long 
era in the early universe where the axion unchains a dynamical screening of the Higgs 
mass. We present a new realization of this idea with the new feature that leaves no 
signs of new physics up to a rather large scale, 10® GeV, except for two very light 
and weakly coupled axion-like states. One of the scalars can be a viable Dark Matter 
candidate. Such a cosmological Higgs-axion interplay could be tested with a number 
of experimental strategies. 



1 Introduction 


Our understanding of Nature is based on the empirical evidence that natural phenomena 
taking place at different energy/distance scales do not influence each other. At present, 
these different phenomena are described by a succession of effective theories with different 
degrees of freedom manifesting themselves as shorter and shorter distances are probed. The 
parameters of the low-energy effective theory are natural if they do not require any special 
tuning of the parameters of the theory at higher energies. 

Wilson [1] and ’t Hooft [2] gave a quantitative meaning to this naturalness principle 
by demanding that all dimensionless parameters controlling the different effective theories 
should be of order unity unless they are associated to the breaking of a symmetry. Numerous 
examples of the naturalness principle to understand the necessity of new phenomena have 
been extensively discussed in the literature (see for instance [3] and references therein). 

The Higgs boson mass and the value of the cosmological constant have been long recog¬ 
nized as two notorious challengers of this naturalness principle, a situation that stimulated 
the creativity of physicists in hnding extensions of the Standard Model at higher energies. 
In most of these efforts to explain the smallness of the Higgs mass, such as supersymmetric 
and composite Higgs models, new physics is predicted to be present at TeV energies. Re¬ 
cently, however, a radically new approach to the Higgs mass hierarchy problem has been 
proposed [4] , in reminiscence of the relaxation mechanism of [5] proposed for explaining dy¬ 
namically the smallness of the cosmological constant (see [6,7] for similar previous ideas). 
In principle, in this new approach no new degrees of freedom around the TeV scale are 
needed anymore to screen the Higgs mass from large quantum corrections. This has of 
course profound implications for the physics agenda of the LHC and beyond. 

Technically, the relaxation mechanism of [4] is based on the cosmological interplay be¬ 
tween the Higgs held h and an axion-like held 0, arising from the following three terms of 
the scalar ehective potential: 

V (0, h) = A^gcf) - + eA^ cos(0//) + • • • , (1) 

where A is the UV cut-oh scale of the model, while Ac < A is the scale at which the periodic 

cos(0//)-term originates and n is a positive integer. The hrst term is needed to force 0 to 
roll-down in time, while the second one corresponds to a Higgs mass-squared term with a 

(positive) dependence on 0 such that diherent values of 0 scan the Higgs mass over a large 

range, including the weak scale. Finally, the third term plays the role of a potential barrier 
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for 0, dependent on /i, necessary to stop the rolling of 0 once electroweak symmetry breaking 
(EWSB) occurs. For this mechanism to work, it is also crucial to have a friction force, coming 
for example from Hubble friction during inflation, in order to make the rolling of (j) very slow 
to access the right minimum during the cosmological evolution. We will discuss the possible 
ultraviolet (UV) origin of Eq. (1) later on. 

At the classical level, the proposed mechanism can be understood in the following way. 
Assuming that 0 starts, at the beginning of the inflationary epoch, at a very large value 
(j) ^ A/g, it will slow-roll until it takes the critical value (pc = A/ at which the Higgs mass- 
squared becomes zero. From this time on, as p continues slowly rolling down, the Higgs mass 
becomes negative, and it is energetically favored to turn on the Higgs held. This raises the 
third term of Eq. (1) up to the point at which p stops rolling. For 1, this occurs for a 
Higgs value v given by 


gA^ ~ 




f 


-e. 


( 2 ) 


This equation arises from demanding that the steepness of the linear (/)-term of the potential, 
hrst term of Eq. (1), equals the steepness of the Higgs barrier, third term of Eq. (1). From 
Eq. (2) we see that we can have v -C A by taking g small enough, which is technically natural 
as g dehnes the spurion that breaks the symmetry p ^ p + 27if. At the quantum level, the 
described cosmological evolution is not much affected provided certain conditions, specihed 
in [4] and discussed later, are fulhlled. Therefore, this mechanism potentially offers a new 
solution to the hierarchy problem. We will refer to this as the cosmological Higgs-axion 
interplay (CHAIN) mechanism. 

For n = 1, the third term of Eq. (1) is linear in h, implying that eA^ must arise from a 
source of EWSB other than the Higgs. This can be the QCD quark-condensate {qq) ~ ^qcd^ 
as proposed in [4]. In this case Ac ~ Aqcd and e ~ ?/„, where is the up-quark Yukawa. 
This model, however, predicts too large a value for the QCD 0-angle, in conflict with neutron 
electric dipole moment constraints. A possible way to £x this problem was explained in [4], 
but it requires a low cut-off scale, A < 30 — 1000 TeV. Alternatively, one could consider 
models in which the condensate comes from a new strongly-coupled sector, a la Technicolor, 
with Ac ~ TeV, or advocate the presence of an additional elementary Higgs doublet. This 
latter case however requires some extra symmetries to keep the second Higgs doublet light, as 
we discuss in Appendix B. These models predict extra physics carrying electroweak charges 
around the TeV scale that can be found in present or near-future experiments. 

For u = 2 on the contrary, the cos{p/ f) term in V{p, h) can arise from the electroweak- 
invariant term cos(0//), where H is the Higgs doublet, and therefore no extra source 
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of EWSB is needed beyond the SM Higgs. As a result, these models can, in principle, allow 
for a larger new-physics scale beyond the SM (BSM). Nevertheless, at the quantum level, 
extra terms can be now induced beyond those shown in Eq. (1). Indeed, just by closing H 
in a loop, we expect, at 0(e), the terms 

eA^cos(0//) , eAlg(pcos{(j)/f), (3) 

to be generated. These terms give a potential to 0 that, unless Ac < n, make it stop slow- 
rolling much before the Higgs turns on. Therefore, if we want the CHAIN mechanism to work, 
we must have again new physics not far away from the weak scale and therefore potentially 
visible in forthcoming experiments. It is important to notice that this new physics is not 
responsible for keeping the Higgs light, but for generating the periodic term of Eq. (1). In 
the particular model of this type discussed in [4], extra fermions were predicted at around 
the weak scale. An important drawback of this type of models is that they must address a 
“coincidence problem”: they must provide a BSM scale Ac that must lie around the weak 
scale with no a priori reason, as the weak scale is determined by Eq. (2). 

The aim of our work is to offer an existence proof that it is indeed possible to devise a 
model that dynamically generates a large mass gap between the Higgs mass and the new 
physics threshold. The proposed model will not have a “coincidence problem” as the only 
new physics scale will be associated with A ~ Ac ^ n. For this to work, we need to make the 
terms of Eq. (3) smaller than the term eA‘^\H\‘^ cos{(j)/f). For this purpose, we will introduce 
another slow-rolling field, a, coupled to cos(0//). During its cosmological evolution, a will 
take a value such that (Tcos(0//) will cancel the terms of Eq. (3). When this occurs, 0 
will be free to move, tracking a downhill. Only when the h-dependent term turns on, 0 
will stop tracking a and reach the minimum hxed by Eq. (2). We will be able to push the 
cut-off scale up to A ~ 10® GeV, providing the hrst example of a natural theory with such a 
large BSM scale. The only new states, 0 and a, will have masses below the weak scale, but 
they will be very weakly coupled to the SM, making them very difficult to detect at present 
and future experiments. Interestingly, as we will see, they could provide the source of dark 
matter needed in the universe. 
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2 Double scanner mechanism 


The key new ingredient of our proposal, with respect to [4], is a second scanning held, that 
we call a. The full potential, up to terms of order e, and g, is given by 

v(,p, + - A" (a - |i/ P + A|i/1* + A(,p, ff, H) cos Wf) , (4) 

where 

A{(j), a, H) = eA^ , (5) 

with 0 < ( 7 , (/cr, e -C 1, while a,j3 and c,p,Ca are 0(1) positive coefficients. We assume that 
all terms of Eq. (4) are generated at the cut-off scale A. For simplicity and clarity, we are 
only considering linear terms in 5 ' 0 /A, but we could have taken a generic function of 5 ' 0 /A 
with the only requirement that it is monotonically decreasing or increasing in a wide region 
of order A/g (and similarly for a with g —)■ g^j). 

From Eqs. (4) and (5), we see that 0 scans the Higgs mass, while a scans A{(j),a,H), 
the overall amplitude -the envelope- of the oscillating term. This dependence of A(0, a, H) 
on a is crucial for our CHAIN mechanism to work, while the other terms in Eq. (5) are 
added since, as we said, they are anyway generated at the quantum level (by loops of H). 
The potential in Eq. (4) is stable under quantum corrections in the small-coupling limit 
( 5 ', 5 'o-,e 1) we consider. A possible UV origin of the periodic term in Eq. (4) is given in 

Appendix A. 

We will study the time evolution of 0, a and H during the inflationary epoch. Inflation is 
needed, as in [4] , to provide the friction that makes the helds slow-roll and reach the desired 
minimum. The time evolution of a is quite trivial, as for e 1, it simply slides down: 

a{t) = cTo - gaAH/{3Hi) . ( 6 ) 

In the cosmological evolution of 0 we can distinguish four stages, depicted in Fig. 1, that we 
qualitatively describe next: 

1) At the start of inflation we assume 0 ^ A/ 5 ' and a > A/g„ such that the Higgs mass- 
squared and the amplitude A are positive. The held 0 is stuck in some deep minimum 
coming from the y4cos(0//) term of Eq. (4), while the Higgs held value is zero. 

11) As a evolves down, the amplitude A decreases until the point at which for 0 the 
steepness of A cos(0//) is smaller than the slope coming from the linear term of Eq. (4), 
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Figure 1: Left; Scalar potential in the {4 >,(t} plane. The hand without harriers is in green 
while the harriers getting high(er) are dark(er) hrown. The blue line shows a possible slow-roll 
cosmological trajectory of the fields during inflation. The dashed purple line is the critical 
line for EWSB. Right; Classical time evolution of (blue curve) in the potential on the 
left. The black lines show the extremal points of the potential, with closely spaced minima 
(bold) and maxima (thin) alternating. (Arbitrary units and scales in both plots.) 


and 0 can start to move down. The region in field-space at which this occurs is shown 
by a “green-band” in Fig. 1. In this region, the bumps from Acos{(j)/f) are very small 
and, for Pa 9, 4> goes down tracking a: 0(t) ~ const. -|- c^g„a{t)/{Cfi)g), which is the 
solution of A r; 0 (this solution neglects effects of size A0 ~ / which correspond to 
the stepwise behavior visible in Fig. 1). 

Ill) When 0 crosses the critical value fie = the Higgs mass-squared term becomes 

negative, turning on H. This gives, according to Eq. (5), a positive contribution to the 
amplitude A, that, for certain values of the parameters of Eq. (4) to be specified later, 
makes the direction of the green-band change as shown in Fig. 1. The field fi cannot 
continue any longer its evolution along the tracking trajectory that would bend in the 
opposite direction in a and would not be compatible with the cosmological evolution 
of a that imposes da/dt < 0. So fi crosses the green-band. 
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Figure 2: Sketch of the four stages in the evolution of (f, marked by the blue dot, in the time- 
dependent effective potential for cf obtained after integrating out a and H but corresponding 
to the same potential as in Fig. 1. (Arbitrary units and scales.) 


IV) Finally, 0 reaches the other side of the green-band (where A becomes large again) and 
gets stuck in a minimum from Acos{(j)/f) as in the model of [4], The held a continues 
going down, making A grow until a hnds its own minimum. 

For the same potential shown in Fig. 1, we show in Fig. 2 four snapshots of how 0 evolves 
in the time-dependent potential 14(0) = 14(0, (j(t), h(0)), obtained after integrating out a 
and h. We are choosing four representative time values corresponding to the four stages 
I-IV. At the stage I and II, this potential has two A 0 regions moving towards each other 
(as indicated by the arrows), that merge at stage III and disappear at stage IV. 

To understand under which conditions the potential of Eq. (4) has the shape shown in 
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Fig. 1, we start by finding the 0 values, called 0*, for which the steepness from Acos{(j)/f) 
is smaller than the steepness from the A^g(j) term of Eq. (4). These are determined by 

< gA^ , (7) 

where we are working in the limit g,ga ^ 1- In this range of 0 values, dV/d(j) = 0 has no 
solutions, leading to a “sliding” region, the minima-free green-band of Fig. 1. 

The value of the neutral Higgs component, (H) = (0, h/\/2), depends on <p according to^ 

- y “ x) ’ ioT(j)<(j)c, (8) 


and zero otherwise. Solving Eq. (7), we obtain the interval(s) 


, 1 C(T 9cF 

<Pc + 

(a - dc) ± — , 

(for 0* > 0c 

C(f) g 

Cfj) € 


1 . C(7 9o- 

0C+ , 

{a - aA ± ^ , 

(for 0* < 0c 

c<f>9 




( 9 ) 


with C 0 = C 0 — 1/(2A). By continuity, both solutions merge at 0* = 0c at a particular value 
o-c = (5'cXc + (3A)/{c„g„)A 

In order for 0(t) to track down a(t), or what is equivalent, for 0(t) to stay in the 0* 
region of Eq. (9) until reaching 0c, the gradient of the dynamical trajectory in the {0, a} 
plane^ inside the 0 * interval, 

d(j){t)/dt _ g 
da{t)/dt ga ’ 

should be larger than the gradient d(f)^:/da of the green-band itself. Otherwise 0(t) crosses 
the green-band too early and gets stuck at some minimum before getting to 0c. From Eq. (9), 
this condition, d(f)(t)/da(t) > dcfy^/da, leads to the requirement 


^4)9 ^ ■ 


( 11 ) 


On the contrary, once 0 < 0c we must demand 0 to exit the green-band (getting trapped 
in some vacuum precisely as needed to explain the smallness of the electroweak scale) which 

^To simplify the discussion, we are ignoring a term — (eA^/A) cos(0//) in the RHS of Eq. (8), that causes 
a transitory period in which the Higgs mass term switches sign repeatedly before stabilizing as negative. 

^In fact, the kink in the sliding-band is located around the point {(j>c^ crA- In terms of these coordinates 
one has A = eA^[c^g{(j) - be) - CagaA - CTc) -b \H\^ / A\. 

^Inside the b* interval, we have A « 0 and b’s slow-roll is driven by the linear term of Eq. (4). 
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implies d(l){t)/da{t) < d(j)^/da and leads to 


■ 


( 12 ) 


This condition is easily satisfied for c'^ < 0, which is equivalent to 2Ac0 < 1, and corresponds 
to the situation depicted in Fig. 1 in which the green-band flips direction at 0 = 0c- For 
2\c^ > 1, the green-band does not switch direction at 0c but its slope changes. This occurs 
however for a small range of where Eq. (12) is also fulfilled: 


2A ^ ^ 2A 


Ca9o 


(13) 


This region of the parameter space can then also provide an explanation to the smallness of 
the electroweak scale. 

Let us emphasize that the CHAIN mechanism described above works independently of 
0j, the initial condition for 0, provided only that 0^ at the initial time lies in the region 
0c < 0i < 0*(ti), which is a natural and sizable range of the available field space. A gfja\H\^ 
term in the potential could spoil the CHAIN mechanism, as the late evolution of a, after 
0 is settled in its minimum, would change the value of the Higgs mass. Therefore we find 
that this term, if present, must be further suppressed by an extra factor e. This does not 
pose any problem to our model as a radiatively generated term arising from Eq. (4) 

by loops of 0, is indeed very small, of 0{e^ga). 


3 Consistency requirements for a small weak scale 

The potential in Eq. (4) involves two scales, A and /, and three small couplings, g,ga- and 
e, apart from a few dimensionless parameters of order unity (including the Higgs quartic 
coupling. A). The cosmological evolution of our model can be broadly described by these 
parameters and two external quantities fixed by the infiaton sector: Hj, the value of the 
Hubble parameter during inflation, and Ne, the number of e-folds. These various parameters 
have to satisfy several constraints in order to provide a natural solution to the hierarchy 
problem, i.e., to dynamically ensure a stable separation between the weak scale and the 
high-energy scale A from which the potential Eq. (4) emerges as a low-energy description. 
These constraints arise from demanding 

1. Dangerous quantum corrections to the potential are kept small. At least at the quantum 
level, additional terms of O(e^) can be generated in the potential of Eq. (4) that could 



potentially spoil our CHAIN mechanism. For instance, terms like^ cos^(0//) or 
e^A^ 5 f 0 cos^( 0 //), depend quadratically on cos(0//), and therefore their amplitudes 
cannot be cancelled by a simultaneously to Acos(0//). These terms are dangerous 
since they give a barrier to 0 at values that can be above the critical 0c- To make sure 
that they remain subdominant to the Higgs barrier of Eq. (4), we must demand 

6 < nVA^ - (14) 


This condition also ensures that the contribution to the Higgs mass coming from the 
eA^lTTp cos(0//) term in the potential is at most of electroweak size and does not spoil 
the tracking behaviour (see footnote 1). 


2. 0 must be trapped by the Higgs barrier. As in the model of [4], the nonzero Higgs held 
must be the only one responsible for stopping 0 from sliding any longer. This is the 
requirement in Eq. (2) that, for our case n = 2 and A = Ac, reads gA^ ~ eA^n^//. 
This can be used to obtain the electroweak scale as a prediction from the model in 
terms of microscopic parameters: 


V 


2 




9 ^ 

e 


(15) 


We will also use this relation later on to get rid of e in terms of the other parameters 
of the model, so that the electroweak scale is reproduced correctly. 


3. Inflation is independent of the 0 and a evolution. We assume for simplicity that 
inhation is driven by another held, the inhaton, that does not receive any back-reaction 
from the evolution of 0 and a. This is possible under the condition that the typical 
energy density carried by 0 and a remains smaller than the inhation scale, i.e.. 


A2 


< Hj 


(16) 


where Mp is the reduced Planck mass, Mp ~ 2.4 x 10^® GeV. 


4. Classical rolling dominates over guantum jumping. We are assuming that the cosmo¬ 
logical evolution of 0 and a is dominated by classical physics. It is therefore essential, 
for the consistency of our solution, that during the cosmological evolution of our system 
the quantum huctuations of the helds, typically of size Hp remain smaller than the 

^See Appendix A for the possible origin of these terms in a particular UV completion. 
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Figure 3: Parameter space for a successful solution of the hierarchy problem ensured by the 
cosmological evolution of the fields cf and a. We have taken A = f and g„ = In 

the upper shaded gray region, the guantum corrections to the potential would not drive the 
system towards the weak scale while in the lower shaded gray region, the two fields would not 
follow their classical paths. Some contour lines of the expansion parameter e (continuous 
black) and of the physical masses of a (dotted blue) and 0 (dashed red) are reported. In the 
shaded red region, the field 0 decays after BBN and has a lifetime shorter than the age of the 
universe. The brown shaded region is excluded by gamma-ray data. Finally, in the shadded 
cyan region on the bottom right corner, the energy density stored in the oscillations of the 
field a exceeds the DM energy density. 


classical field displacements over one Hubble time, i.e., Aa Hj ^da/dt ~ Hj "^dV/da 
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for the case of a. This condition for classical rolling [4] simply reads 


Hf 




( 17 ) 


for a. Due to Eq. (11), the classical-rolling condition for 0 is automatically guaranteed 
whenever Eq. (17) is fulhlled (we assumed that c<^ ~ Co- ~ 1). 


5. Inflation lasts long enough for complete scanning. In order for the above mecha¬ 
nism to work for generic initial held conhgurations, it is essential that the range 
scanned by 0 and a during inhation be of order or larger than A/g and A/g„ re¬ 
spectively. This is ensured by requiring a long enough period of inhation, namely, 
NgAa ~ NeHj'^dV/da > A/g^, which leads to 


N >AL 

e ~ _9 A 2 




( 18 ) 


The typical duration of the diherent stages of the cosmological evolution of 0 is of the 
same order, with the exception of stage III, which is much shorter, of order ~ 
(^/M)A^e- 


Combining these various consistency conditions, together with Eq. (11), we obtain that the 
couplings g^ and g are bounded to the interval: 

^ 9(7 ^ 9 ^ /A^' 

Since / cannot be much smaller than A, as this latter is the scale at which the cos(0//) 
term is generated, we obtain from Eq. (19) an upper bound on the cut-off of our model; 

A < (n^M|)^/^ ~ 2 X 10^ GeV. (20) 

The bound in Eq. (19) dehnes the region of the parameter space of the model consistent 
with a natural solution to n A. This is shown in Fig. 3, where, for concreteness, we have 
taken / = A and ga^ = O.l^f. 


4 Quantum spreading 

The discussion of the cosmological evolution in Section 2 was fully classical. As we saw 
in more detail in the previous Section, the model parameters can be chosen to ensure that 
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the slow-roll of the helds 0 and a is indeed dominated by classical evolntion, with qnantnm 
flnctnations (driven by Hj) playing a snbdominant role. As discnssed in [4], however, one 
shonld pay attention to the snbleading qnantnm spreading to ensnre that it does not spoil 
the solntion to the hierarchy problem. 

This issne can be attacked nsing a Langevin eqnation (that adds to the nsnal classical 
evolntion of the held a noise term that captnres the random hnctnations controlled by Hj) 
or nsing the Fokker-Planck eqnation® 

= T ^ (Hip) (21) 

dNe d(j) V3-^/ d(j)) d(j)‘^ ySvr^ / 


for Nf.) the probability (density) of hnding the held in the interval (0, (j) + dcp) in a 
given Hnbble patch after e-folds of inhation. The hrst term in the eqnation describes the 
ehect of classical slow roll while the second describes the qnantnm hnctnations indnced by 
inhation. 

It is convenient to dehne the average valne of 0, (0), and its variance = {{4>— (0))^), 
where 

/ CO 

P{ct>.N:)F{ct>)dct> . ( 22 ) 

■CO 

Using the Fokker-Planck Eq. (21), one gets 


dN, 


1 

“3^ 

47|-2 


dV\ 

Ml ’ 


3H] 





(23) 

(24) 


The hrst eqnation shows that (0) slow-rolls classically down an averaged potential slope. The 
second eqnation describes the qnantnm spread of the held. Besides the positive contribntion 
that makes the variance grow with if/, as A0 ~ there is an ehect from the potential 

that can eventnally stop the qnantnm spreading {e.g. in a potential well). 

We can nse these eqnations to gain some qnalitative nnderstanding of how the held 0 will 
spread dnring its cosmological evolntion (described in Section 2). Dnring stage I, 0 is stnck 
in a deep minimnm where we can neglect qnantnm spreading ehects. In stage II, 0 rolls at 
the bottom of the 0* range and it is a good approximation to assnme that (0) ~ (0*)min, 
the low end of the 0* region. The potential aronnd that point is asymmetric, with a nearly 
constant and very small slope in the 0* range and a steeper slope ontside that range (the 

®For more details on the use of these equations to describe fluctuating light fields (like the Higgs) during 
inflation, see e.g. [8]. 
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hilly side that actually slows down the held so that it tracks (0*)min)- As a result, the spread 
of 0 is also asymmetric, being suppressed in the hilly side and larger inside the 0* range. 
The latter spread starts growing at the beginning of stage II and, depending on how large 
Hi is and how long stage II lasts, it might reach a steady state (with dAcj)"^/dNf> ~ 0). 
Irrespective of that, it is guaranteed that the hnal quantum spread will be at most of the 
same order of the size of the 0* region (A0* = //(c<^e) ~ //e) because the large barriers 
beyond that range suppress the spreading.® How far away from the 0*-region the barriers 
will stop the spreading can be estimated (by comparing the classical falling speed with the 
uphill quantum jumps) as 

(25) 

which is the width of the band times a factor smaller than 1, see Eq. (17). The same applies 
to stage III, with a similar bound on the hnal spread A0 < //e. In terms of the Higgs 
held expectation value, using Eq. (8), this translates into an spread of the order g /A/(2Ae), 
which is of the order of as we saw in the previous Section, Eq. (15). So the mechanism to 
solve the hierarchy problem is not spoiled. As in [4], our modihed CHAIN mechanism does 
not oher a solution to the cosmological constant problem. 

Quantum spreading will also ahect the held a, which rolls down a much simpler potential. 
Its spread will grow as Aa ~ y/N^Hj and its hnal value at the end of inhation will depend 
on the shape of the potential near the hnal a minimum. Concerning the evolution of the 
held 0, the spreading of a will not matter (as long as classical rolling dominates): 0 will 
simply track a wherever it is in a given Hubble patch. 

5 Cosmological signatures 

As we saw in the previous Sections, in contrast with the models of Ref. [4] , the new-physics 
scale of our model can be as large as A ~ 10® GeV, and therefore we do not expect any new 
state around the weak scale. Only the two additional scalars a and 0 are lighter than, or at 
most around, the weak scale. However, as we will see below, they are very weakly-coupled to 

®There is always a small probability that the field leaks to neighbouring minima (with the overall slope 
biasing the leakage towards lower minima) but the exponential suppression (^ exp[—87r^y/(3i4j)]) of the 
process overcomes the large number of e-folds in stage IV. Notice also that the minimum at which 0 ends up 
is extremely stable against decay by quantum tunneling after inflation is over. The tunneling bounce action 
can be estimated (in the thin-wall approximation) to be S 4 ^ e^f/{Ag^) > (A/u)®, which leads to a vacuum 
lifetime many orders of magnitude larger than the age of the universe. 
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the SM states in most of the parameter space, and thus can only have some phenomenological 
impact through astrophysical and cosmological effects. 


5.1 Properties of 0 and cr 


We start by deriving the properties of the 0 and a scalars. After the slow-rolling process 
ends and a settles in a minimum, no cancellation is expected in the ^(0, a, H) amplitude, so 
that A{(j), a, H) ~ eA^. The mass of (j) is thus controlled by Acos(0//) and can be estimated 

A5 

(26) 


2 


P 


A' < 2 


where we used Eq. (14) and Eq. (15) to obtain the second equality and the upper bound on 
m^. For (j we expect that higher-order terms in g^a/K, not shown for simplicity in Eq. (4), 
give it a mass of order 

ml ~ glk^ < ml . (27) 


In the allowed part of the parameter space of our model the masses of the two scalars 
can change by many orders of magnitude, spanning the range [10“^°, 100] GeV for 0 and 
[10“^^, 10“^] GeV for a. Gontours of constant m^ and rrio- are shown in Fig. 3. 

These two scalars interact with the SM particles mainly through a mass mixing with the 
Higgs. The corresponding mixing angles can be estimated as 
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(28) 


Notice that the (p — h mass mixing coming from ~ eA^(n//) sin(0//) is suppressed at 
the minimum where we have sin(((;/))//) 9flW ~ !iP lA The first contribution 

in 6(jh arises at tree-level, whereas the second one originates from a 0-loop. For most of the 
parameter space we consider, this loop term dominates over the tree level one. The scalar 
potential Eq. (4) also gives rise to interactions between 0 and the Higgs, not suppressed by 
the small mixing angle 9^h, that are of order 


(pfphh : eA^/p , (jxph : evA^/p , (29) 

^This is to be contrasted with the beginning of Phase IV when sin(0//) ~ 1, as used to derive Eq. (15), 
since barriers are smaller at this earlier stage. At the end of Phase IV the barriers have grown large, and 0 
is close to the minimum of its cosine potential. 
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and will play an important role in the thermal prodnction of 0. The decays of 0 and a are 
mediated by the mixing with the Higgs, and thns the widths are given by 

^^ h{m^) , T^ ^ elfj:h{ma), (30) 

where Th{mi) is the SM Higgs width evalnated at rrih = mi. Contonrs for T^o- are shown 
in Fig. 3 (the valnes of the width Thimi) are snbject to large theoretical nncertainties in 
the mass region aronnd 1 GeV where several hadronic decay channels open np [9]; we nsed 
the expressions given in Ref. [10] -see also Refs. [11,12]). For masses below 2me ~ 1 MeV, 
we have Thimi) ~ (mi/mh)^ Th^-yy{mh), and therefore, in a major part of the parameter 
space, 0 and a have snppressed decay widths controlled by the decay into photon pairs. As 
shown in Fig. 3, there is a sizable part of the parameter space in which 0 is cosmologically 
nnstable (F,^ > Hq, where Hq is the present Hnbble valne), bnt snfficiently long-lived to 
decay after Big Bang Nncleosynthesis (BBN) (F^ < Hbbn = hf(T = 1 MeV)). As we will 
see in the following, this region of the parameter space can be constrained by cosmology. On 
the other hand, a is cosmologically stable in most of the relevant parameter space, and can 
decay within the age of the nniverse only in a small corner of the parameter space, namely 
for Qcr > 10“® and A < 10^ GeV. 

We can now easily estimate the cosmological abnndances of 0 and cr, either stored in 
late classical oscillations (vacnnm misalignment) or from thermal prodnction. This will 
allow ns to set bonnds on the model from overclosnre of the nniverse, post-BBN decays or 
astrophysical constraints. 

5.2 Impact of (f) and cr on standard cosmological predictions 

In this work we assnme for simplicity that, once both 0 and a have settled in their minima, 
inflation ends with an nnspecihed reheating period. We will assnme a reheating temperatnre 
higher than the EW scale in what follows. 

Abundances of <p and cr from vacuum misalignment 

If after inflation and reheating, the helds 0 and a end np displaced from their minima, 
they will fall towards them, oscillating aronnd them if their lifetimes are large. The energy 
density stored in the held oscillations behaves like cold dark matter and can potentially 
overdose the nniverse today or dissociate light elements if the decay takes place dnring or 
after BBN. At the start, the held energy density is dominated by the potential energy, bnt 


15 


as the fields roll to their minima and pick up speed, the kinetic energy grows. When both 
contributions are of similar size, one can properly talk about oscillating fields. This happens 
for the Hubble value if < mj (i = 0, a) or, equivalently, in a radiation-dominated universe, 
at ~ \/miMp. In that oscillating regime, the energy density scales as that of non- 
relativistic matter, Pi(T) ~ plniiT , where is the initial amount of energy at the 
start of the oscillating regime (this is smaller than the original potential energy density, but 
only by an order one factor). 

Let us start considering the a field. We expect that during inflation a slowly rolled 
down to its global minimum, somewhere in its ^ K/ range, as this requires a number of 
e-folds similar to the W estimated in Eq. (18). Because of quantum effects (see Section 4), a 
reached the minimum with a spread y/N^Hj. We can use this result to estimate the typical 
displacement from the minimum at the end of inflation: 

{Aa)ini ~ VNeHj . 

This displacement is quite large though still smaller than A/g^j, and hence we can estimate 
the amount of energy density as ~ Using ~ Hj/{A^gl), we finally get 

~ Hj. The energy density stored in a oscillations today, relative to the critical energy 
density, is then, using Eq. (16), Ho- ~ (4 x (A/10® GeV)^^^^ . The bound to avoid 

universe overclosure (Ho- < 1), translates into 

/ A \ 

The contour H^ = fipM is plotted in Fig. 3, limiting the excluded blue region. It shows that 
a can be a good dark matter candidate in the region where the bound (31) is saturated, in 
particular at large A. For certain values of m^, there can be other cosmological constraints. 
For example, for Ho > H£)m/ 20, the mass range 10“®^ eV < < 10“^^'® eV is excluded by 

structure formation [20], while masses around m„ ~ 10“^^ eV may be constrained by Black 
Hole superradiance [14]. Interestingly, for the particular case ~ 10“^^ eV, a could be 
searched for by the SKA pulsar timing array experiment [19]. Let us finally notice that there 
are ways to go around the bound (31), for instance, by assuming a late entropy production 
after a has started to oscillate, as can occur if reheating is a very slow process such that 
Trh < T^sc [ 13 ]- 

The situation for 0 is rather different. At the end of its evolution, cj) is trapped in 
a region with high barriers and its displacement from the minimum originates from the 
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possible quantum spreading. The initial energy density arising from this displacement was 
at most ~ Hf, that, since 3> rUo- and then gives today a completely 

negligible effect. 


Thermal production of <p 


Thermal production of 0 arises mainly from the couplings of Eq. (29). In particular, from the 
00/i/i-couphng we can have double-production from the thermal bath via h + h ^ ® At 

T > ruh, this double-production cross-section is estimated to be (uav) ~ e2(AV/4)/T2. This 
implies that 0 can reach thermal equilibrium only for T in the interval [rrih, e^Mp{A/f)^], in 
which the 0 production rate is faster than the rate of expansion. This region corresponds 
roughly to the area above the T^ = Hbbn line of Fig. 3, so we conclude that in most of the 
parameter space, 0 never thermalizes.® 

The number density of 0 produced thermally is obtained by solving the Boltzmann equa¬ 
tion 

^ + 3Hn^ = - (aAv) (nj - , (32) 

where Utp^eq is the equilibrium number-density of 0. This equation can be conveniently re¬ 
written in terms of the dimensionless quantities x = m^/T and = n^/s, where s is the 
entropy per comoving volume, s = 27r^p*sT^/45. Assuming a radiation-dominated era, with 
energy density pp = 7r^p*T‘^/30 (here, p* ~ ~ 100 counts the number of relativistic 

degrees of freedom) and using that 10, ^<t>,eq i'^ large portion of parameter space in 

which 0 does not thermalize, one gets: 


dl0, ^ {aA^)Cm^Mp y2 
dx x‘^ 


(33) 


where C = 27rv^fi'*s/(45y^) ~ 13.7. For relativistic 0, a; 1, the equilibrium density is 
approximately given by 10.q ~ 0.278/p*5. This leads to the approximate formula 

U(T) ~ . (34) 

^Double production can also be mediated by the process SM + SM -A h^*'> -A cj) + cj) induced by the 
00/i-coupling, which can lead to a similar thermal production as the one discussed here. Single production, 
on the other hand, is due to interactions that are linear in the 0 field and are thus suppressed by the small 
mixing angle and can be neglected. 

^This also implies that we can neglect thermal corrections to the potential for 0 in the analysis of its 
cosmological evolution. 
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The 0 production is maximal at T ~ ruh- Using Eq. (34) evaluated at T ~ mh and assuming 
no late entropy production, we can deduce, in the parameter region where 0 is cosmologically 
stable, the contribution of 0 to dark matter today, ~ m^Y^so/pc where Sq is the present 
entropy density. We find that T^ is subdominant; for example, along the contour = Hq 
in Fig. 3, we have <3 x 10“'^, while for T,^ ~ we hnd < 7 x 

Constraints from BBN and Gamma-Ray observations 

As mentioned earlier, there is a region of parameter space in which 0 is not cosmologically 
stable and decays after BBN. This is problematic if the decay of 0 injects into the thermal 
bath an energy per baryon Ep,b > O(MeV), leading to a modihcation of the predictions for 
the abundances of the light elements. Since Ep b ^ where is the abundance 

prior to decay (i.e. the value of Eq. (34) obtained as if 0 were stable), this results in the 
bound < 10“^^ GeV, which however could be weakened sensitively depending on the 

precise value of the lifetime [15]. In the parameter space region where 0 decays between 
cosmic times ~ 1 s and ~ 10^”^ s, varies from ~ 100 GeV down to 1 MeV and 
varies from 10“^ to 10“^^ GeV. Above the di-pion threshold, hadronic decays constrain short 
lifetimes T^^ < 10^ s while for lighter 0, electromagnetic decays constrain large lifetimes 
> 10^ s [16,17]. In addition, the Gosmic Microwave Background (GMB) constrains lifetimes 
~ [10^'^ — 10^^] s for Ep,b down to 0(eV). Therefore, we expect that most of the region of the 
parameter space delimited by the lines = Hbbn and T,^ = Rq ia Fig. 3 is excluded. A 
dedicated analysis is however needed to derive the precise excluded regions that is beyond 
the scope of this paper. 

On the other hand, for regions in which the 0 lifetime is larger than the age of the 
universe, there are strong constraints coming from decays generating a distortion in the 
galactic and extra-galactic diffuse X-ray or gamma-ray background. In particular, sub-GeV 
dark matter decaying into photons should satisfy tdm ^ 10^’^ s [18]. Since the gamma-ray 
flux scales as d^-y/dE oc V^T^, we can translate this bound into > 10^'^ s x Q^/Qdm- 
This excludes the thin brown band of Fig. 3. 

We stress again that the cosmological constraints derived above can be evaded if the 
temperature of the universe never reaches mh, in which case the thermal production of 0 is 
suppressed. 
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6 Conclusion 


The aim of this paper has been to provide an existence proof of a model, based on the idea 
of [4], that can naturally accommodate a small electroweak scale without (sadly, to say) 
requiring visible new-physics at present and far future colliders. The model is based on a 
cosmological evolution of the Higgs and two axion-like states whose back-reactions lead to 
a naturally small electroweak scale. The only new-physics of the model consists of these 
two scalars, (j) and a, that in most of the parameter space are very light and very weakly 
coupled to the SM. Therefore strategies to detect them are completely different, as they do 
not require powerful high-energy colliders, but dedicated searches in the sub-GeV regime. 

Interestingly, one of the two light states, namely a, could be a dark matter candidate. 
Because of its small mass, a starts oscillating very late and can therefore make a large 
contribution to dark matter today. The field 0 cannot contribute to more than 
For this maximum value, it may still be detected in gamma-ray observations from its late 
decay. 

Part of the parameter space of our model can be tested through observations of the 
diffuse gamma-ray backgrounds, black hole superradiance and even in pulsar timing arrays. 
In addition, there is a rather rich BBN and CMB phenomenology which motivates a more 
thorough study. Unfortunately, hfth-force signals and Equivalence Principle violations in 
intermediate mass ranges seem too small to be seen in the near future. 

There are important discussions that we leave for the future. An interesting aspect of 
our mechanism is that it does not depend on initial conditions for the new scalars, provided 
that it takes place during a long period of inflation. Similarly to [4], we need a very long 
period of inflation, > 10^®, and an inflation scale rather low, < 10® GeV. In addition, 
inflation requires super-Planckian held excursions. All of these issues deserve further study. 
It seems clear that the possibility to further increase the effective cutoff A might help with 
these questions. 

To close on a more philosophical note, the ideas proposed in [4] and pursued here represent 
a new twist in the long and fruitful history of the interplay between particle physics and 
cosmology. While in the past particle physics has been a crucial ingredient to understand the 
cosmological history of our universe, if these new ideas were correct, cosmological evolution 
would be a key ingredient in the understanding of some key parameters of particle physics. 


19 


Acknowledgments 


We thank Joan Elias-Miro, Ednard Masso and Marc Riemban for very usefnl discussions 
in the early stages of this project and David E. Kaplan, Kai Schmidt-Hober and Alexan¬ 
der Westphal for very interesting discussions. We acknowledge support by the Spanish 
Ministry MEC under grants FPA2014-55613-P, FPA2013-44773-P, and FPA2011-25948, by 
the Generalitat de Catalunya grant 2014-SGR-1450 and by the Severo Ochoa excellence pro¬ 
gram of MINECO (grant SO-2012-0234). CG is supported by the European Commission 
through the Marie Curie Career Integration Grant 631962. CG and GS are supported by 
the Helmholtz Association. The work of AP has been supported by the Catalan ICREA 
Academia Program. 

A A partial UV completion of the model 

Let us start the discussion by summarizing the basic building blocks of the potential of 
Eq. (4). The hrst ingredient is the effective potential for the scalar 0 that contains two types 
of terms. The periodic term,^° associated to a discrete shift invariance, 0 —)■ 0-7 27r/, and the 
linear terms that break this symmetry by the small coupling g. It is crucial that both types 
of terms couple to for the CHAIN mechanism to work. Due to Eq. (15), we have g e 
and therefore the linear terms are always smaller than the periodic term. In addition, the 
model also contains a second held, a, whose interactions are controlled by the small coupling 
ga, in such a way that a always appears in the combination g„a/A. In analogy to the 0 held, 
we can associate a to a continuous shift invariance, whose breaking is controlled by g^. It is 
crucial to have a coupling of a to the periodic function of 0, but not a direct coupling to the 
Higgs that would spoil the CHAIN mechanism, as we mentioned before. Quantum ehects 
via a loop of 0-held can of course generate this coupling, cr|i7p, but its size of order 0{e‘^g„) 
is small enough and does not destabilize the weak scale. 

Although each of these terms needs an appropriate UV origin, here we are only interested 
in providing a UV model that explains the generation of the periodic terms in the potential 
of Eq. (4). The model is a generalization of the non-QCD model proposed in Section 3 of 
Ref. [4]. It consists of a gauge sector, which, analogously to QCD, we take to be SU{N), 
with two Dirac fermions, L and N, belonging to the fundamental representation of SU{N). 

^°The form of the periodic function is not important for the model. The choice cos(0//) has been selected 
in view of a possible axion-like UV completion, as we will see in the following. 
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N N 

Figure 4: Left: Diagram generating (f)NN at the radiative level. Middle: Diagram con¬ 
tributing to the coupling NN\H\'^. Right: Diagram generating an O(e^) contribution to 
(NNf. 



Under the SU{2)l x U{1)y SM group, L has the quantum numbers of a lepton doublet, while 
is a singlet. We assume that the SU{N) gauge sector becomes strongly-coupled at the 
scale A. A key ingredient of the model is the presence of a specific set of mass and interaction 
terms for the fermions that break the accidental global symmetries. We assume that the L 
and N fields have Dirac masses (here and in the following we neglect 0(1) parameters): 

h^raass — d^LL cANN , (35) 

and couplings to the SM Higgs given by 

Cvuk = Ve LHN + h.c.. (36) 

Finally, interaction terms of the singlet N to the a and 0 fields are included with couplings 
of order eg and eg^j respectively 

Cn = egcpNN -|- eg^crNN . (37) 

As can be seen from the Lagrangian above, we have associated to each N field a coupling 
^/e -C 1. In the limit e —)■ 0 the theory acquires an additional chiral invariance (broken only 
by the axial anomaly). It is interesting to notice that even if we do not introduce in the 
Lagrangian the coupling of the 0 field to N, it is nevertheless generated at the radiative level 
due to the presence of the gA(j)\H\‘^ coupling in the effective Lagrangian, as shown by the 
left diagram of Fig. 4. 

We also assume that the 0 field interactions are invariant under a shift-symmetry, 0 —)■ 
0 -|- c, up to the explicit breakings due to g, and an anomalous interaction term 

, ( 38 ) 
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where G' denotes the SU{N) field strength. Analogously to the QCD axion, the field 0 
acquires a periodic effective potential as a consequence of the chiral anomaly. The best way 
to estimate the size of this contribution is to perform a chiral rotation for N such that one 
eliminates the term (38) but generates rn^NN —)■ where uin is the effective 

mass of N: 

A + g^a + g4> -, (39) 

with the last term arising from the middle diagram of Fig. 4. Using (NN) ~ A^, the term 
NN + h.c. gives, at 0(e), 

U ~ A^rriAr cos(0//). (40) 

Equation (40), together with Eq. (39), reproduces the periodic term of Eq. (4). 

Using this explicit UV model we can also analyze possible additional contributions to the 
0 potential. For example, at O(e^) we expect contributions to the 0 potential coming from 
terms {NNY generated at the quantum level from diagrams as the one shown in Fig. 4. After 
the chiral rotation described above, that eliminates Eq. (38), we have (NNY —)■ 
which leads to a term for the 0 potential ~ e^A^cos^(0//). As we discussed in the main 
text, in order for the relaxation mechanism to work we need to suppress these terms with 
respect to the leading potential in Eq. (4). This leads to the condition in Eq. (14). 

B Relaxation in a two Higgs doublet scenario 

In this appendix we present a relaxation model based on a two Higgs doublet model (2HDM). 
The motivation for this is, as mentioned in the introduction, to generate the term /icos(0//) 
that requires a second source of EWSB. If the second Higgs is also elementary, we must find 
a way to keep its mass also small. Otherwise, we expect that, at the end of the relaxation 
process, only one Higgs being light, while the second doublet having generically a mass of 
the order A. 

To solve this problem we can advocate an additional global SU{2)r invariance at the 
scale A under which the two Higgses transform as a doublet, {Hi,H 2 ), ensuring that both 
have the same masses and quartic couplings. This symmetry guarantees that the masses of 
the two Higgses are canceled simultaneously by the 0 field. The S'U(2)/j symmetry can be 
easily extended to the third quark generation sector by considering a type-H 2HDM in which 
the bji and components form an S'U(2)j:j-doublet, and the Yukawa term is given by 

yQL{H^ H 2 ){bRtRf + h.c.. (41) 
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Below A, the most important source of breaking of the SU{2)ji invariance arises from the 
SM hypercharge gauge-boson that of course distinguishes from (and Hi from H 2 ). 
The effective potential involving the Higgses and the held 0 reads 

+ km^\Hi\^ + + h.c) + • • • , (42) 

where m|^(0) = —A^ (1 — g(f)/A). In the above expression we included the leading SU{2)ji- 
breaking that arises from two-loop diagrams involving fermions and the hyper¬ 

charge gauge-boson, and whose size can be estimated to be 

( 43 ) 

If we want this model to have two independent sources of EWSB, Hi and H 2 , we encounter 
a phenomenological problem. The model will contain extra light Higgs, of mass of order 
< '/Xv, that have not been seen at the LHC. For this reason, we are forced to look for 
2HDMs in which only one Higgs, H 2 , is responsible for triggering EWSB, while the other 
Higgs, Hi, gets a vacuum expectation value from its mixing with H 2 . The masses of the 
extra Higgs will then be of order Am that can be heavier than the EW scale and escape 
present detection. 

The dynamics of this model is quite straightforward to understand. When, as a conse¬ 
quence of the slow-roll of 0, becomes negative, only H 2 turns on, as the Hi mass is 

still positive if, as we assume from now on, m|^(0) -|- Am^ > 0. Nevertheless, the periodic 
term in the potential automatically induces a tadpole for Hi which displaces it from the 
origin. In such a scenario we have 

{H 2 ) ~ V , {Hi) ~ . (44) 

As the 3rd family Yukawas are equal, Eq. (41), we have mb/mt — {Hi)/{H 2 ), that requires 

eA^/Am^ mb/mt- (45) 

In order for the relaxation mechanism to work, we need to ensure that the potential is 
radiatively stable. The fact that we have two Higgses helps to protect the periodic term of 

additional breaking is induced at one loop by light fermion Yukawas. These effects are however 
subleading with respect to the ones induced by the hypercharge gauging. 
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Eq. (42) that can only be self-renormalized. Therefore the only one-loop extra term that we 
can have is given by 

^^e2A^cos2(0//), (46) 

that should be smaller than the periodic term in Eq. (42) (recall that here, contrary to our 
main model, we are not introducing the held a). This leads to the condition 


e < IGtt 


2mb 
rrit Is? ’ 


( 47 ) 


which, together with Eq. (43) and Eq. (45) gives an upper bound on the cut-off A: 


A < ^ 106 GeV. (48) 

gy 

This model is thus characterized by three scales: the weak scale v, the mass of the second 
Higgs doublet tuhi ~ ~ scale of new dynamics A < 10^ v. This implies 

that only the second Higgs doublet would be reachable at the next LHC runs. 


References 

[1] K. G. Wilson and J. B. Kogut, Phys. Kept. 12 (1974) 75. 

[2] G. ’t Hooft, NATO Sci. Ser. B 59 (1980) 135. 

[3] M. Dine, [hep-ph/1501.01035]. 

[4] P. W. Graham, D. E. Kaplan and S. Rajendran, [hep-ph/1504.07551]. 

[5] L. F. Abbott, Phys. Lett. B 150 (1985) 427. 

[6] G. Dvali and A. Vilenkin, Phys. Rev. D 70 , 063501 (2004) [hep-th/0304043]. 

[7] G. Dvali, Phys. Rev. D 74, 025018 (2006) [hep-th/0410286]. 

[8] J. R. Espinosa, G. F. Giudice and A. Riotto, JCAP 0805 (2008) 002 [hep-ph/0710.2484]; 
J. R. Espinosa et ah, [hep-ph/1505.04825]. 

[9] J. R. Ellis, M. K. Gaillard and D. V. Nanopoulos, Nucl. Phys. B 106 (1976) 292. 

[10] F. Bezrukov and D. Gorbunov, JHEP 1005 (2010) 010 [hep-ph/0912.0390]. 


24 



[11] K. Schmidt-Hoberg, F. Staub and M. W. Winkler, Phys. Lett. B 727 (2013) 506 [hep- 
ph/1310.6752]. 

[12] J. D. Clarke, R. Foot and R. R. Volkas, JHEP 1402 (2014) 123 [hep-ph/1310.8042]. 

[13] G. F. Giudice, E. W. Kolb and A. Riotto, Phys. Rev. D 64, 023508 (2001) [hep- 
ph/0005123]. 

[14] A. Arvanitaki, M. Baryakhtar and X. Huang, Phys. Rev. D 91, no. 8, 084011 (2015) 
[hep-ph/1411.2263]. 

[15] R. H. Cyburt, J. R. Ellis, B. D. Fields and K. A. Olive, Phys. Rev. D 67 , 103521 (2003) 
[astro-ph/0211258]. 

[16] H. An, M. Pospelov, J. Pradler and A. Ritz, Phys. Lett. B 747 , 331 (2015) [hep- 
ph/1412.8378]. 

[17] M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 71 , 083502 (2005) [astro- 
ph/0408426]. 

[18] R. Essig, E. Kuflik, S. D. McDermott, T. Volansky and K. M. Zurek, JHEP 1311 (2013) 
193 [hep-ph/1309.4091]. 

[19] A. Khmelnitsky and V. Rubakov, JCAP 1402, 019 (2014) [astro-ph.CO/1309.5888]. 

[20] R. Hlozek, D. Grin, D. J. E. Marsh and P. G. Ferreira, Phys. Rev. D 91 (2015) 10, 
103512 [astro-ph.CO/1410.2896]. 


25 



